Minimal instances for toric code ground states 



o 

(N 

On! 



Nicolai Lang* and Hans Peter Biichler 
Institute for Theoretical Physics III, Universitat Stuttgart, 70550 Stuttgart, Germany 

(Dated: July 2, 2012) 

A decade ago Kitaev's toric code model established the new paradigm of topological quantum 
computation. Due to remarkable theoretical and experimental progress, the quantum simulation 
of such complex many-body systems is now within the realms of possibility. Here we consider the 
question, to which extent the ground states of small toric code systems differ from LU-equivalent 
graph states. We argue that simplistic (though experimentally attractive) setups obliterate the 
differences between the toric code and equivalent graph states; hence we search for the smallest 
setups on the square- and triangular lattice, such that the quasi- locality of the toric code hamiltonian 
becomes a distinctive feature. To this end, a purely geometric procedure to transform a given 
toric code setup into an LC-equivalent graph state is derived. In combination with an algorithmic 
computation of LC-equivalent graph states, we find the smallest non-trivial setup on the square 
lattice to contain 5 plaquettes and 16 qubits; on the triangular lattice the number of plaquettes and 
qubits is reduced to 4 and 9, respectively. 
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The quest for topological quantum computation [1] re- 
lies on well controllable quantum systems exhibiting any- 
onic statistics. Several theoretical studies have proposed 
microscopic models featuring such topological phases; the 
most prominent example being the toric code [2]. While 
it is generally accepted that the fractional quantum Hall 
effect is well described as a topological state of matter [3- 
5], alternative setups such as cold atomic gases [6] and 
ion traps [7] are considered as suitable building blocks for 
a quantum simulation of topological matter. However, 
such systems naturally appear with a small number of 
involved particles. Regarding actual implementations it 
is therefore a fundamental question, how many particles 
the simulated system must involve in order to observe 
the characteristic properties of the system in question. 
In this letter, we address this question for the toric code 
and present setups with a minimal number of qubits in 
order to obtain non-trivial ground states. 

In the last years great experimental progress led to 
the generation of highly entangled multi-qubit states in 
realistic settings, such as cold atomic gases in optical lat- 
tices [8, 9] , Rydberg states [10-12] , trapped ions [7, 13- 
15] and polarized photons [16-18]. On the theoretical 
side, a powerful tool in characterizing and describing 
states with mutilpartitc entanglement is given by the 
concept of graph states [19, 20], and several realistic se- 
tups for the generation of graph states have been pro- 
posed [21-24]. 

In contrast to the coherent and measurement based 
state preparation above, the concept of quantum simu- 
lation [25, 26] of complex many-body Hamiltonians pro- 
vides an alternative approach for the generation of topo- 
logical phases as a suitable resource for topological quan- 
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Figure 1. (Color online) A: The (non-degenerate) ground 
state of a single toric code plaquette is a GHZ state and equiv- 
alent to the shown graph state. Qubits (o) are located on 
the edges. Vertices ( • ) denote star operators A 3 and grey 
faces plaquette operators B p acting on the adjacent qubits. 
The LC-equivalent graph state connects the qubits directly 
{——)■ B: The 8 neighbours ( ) of the centered qubit 

( • ) in the toric code on a square lattice. Neighbours and 
the centered qubit have a star and/or plaquette operator in 
common. 



turn computation; the paradigm of which has been es- 
tablished by Kitaev's toric code [2]. Subsequently sev- 
eral setups for the quantum simulation of the latter were 
discussed [6, 27-29]. Special focus was laid on the ob- 
servation of mutual anyonic statistics for the elemen- 
tary excitations [30-33] , and experimental success within 
the realm of photonic quantum simulation was reported 
[34, 35]. 

However, the systems considered so far are small and 
their ground states tend to be well known multi-qubit 
states, e.g., the ground state of a single toric code pla- 
quette is a simple GHZ state [36], as shown in Fig. 1 (A). 
Hence we argue that a proper characterization of toric 
code systems implies the comparison of its ground states 
with other known multi-qubit states, such as graph 
states. 

In this letter, we address the question on the minimal 
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Figure 2. (Color online) The minimal non-local toric code setup on a square lattice (A) and on a triangular lattice (B) with 
LC-equivalent graph states. The system depicted in (C) was used to prove the non- locality of the square lattice setup in (A). 
Notation: Qubits ( O ) are located on the edges. Vertices ( • ) denote star operators A B and grey faces plaquette operators B p 
acting on the adjacent qubits. The drawn LC-equivalent graph states connect qubits directly via local ( — ) and non-local 
(••••) edges. The graph states are based on the black spanning trees ( — ), see Definition 1 and Theorem 1 for details. 



system required to observe such characteristic properties 
of the toric code. Let us explain this in more detail: 
The implementation of the toric code naturally gives rise 
to the concept of quasilocal [37] interactions and nearest 
neighbour sites: The qubits are located on the links of a 
graph which is embedded into a surface, see Fig. 1 (B). 
The qubits connected to the same vertex s give rise to 
the quasilocal stabilizer A s = Y\ iGs erf, while the qubits 
surrounding a plaquette p are involved in the many-body 
stabilizer B p = Jl iep erf; the latter are qubits connected 
to a common vertex of the dual graph. This leads to a 
natural definition of nearest-neighbour (or vicinal) qubits 
in the toric code, as qubits involved in the fundamental 
stabilizers A s or B p . For a square lattice, the number 
of nearest neighbours is therefore eight, see Fig. 1 (B). 
It is a fundamental property of the toric code, that it is 
generated by quasilocal stabilizer operators, and there- 
fore, its experimental realization requires only quasilocal 
interactions. 

As mentioned above, for small systems the ground 
state of the toric code is equivalent to well known states: 
for example, the non-degenerate ground state of the setup 
with four qubits depicted in Fig. 1 (A) is equivalent to 
a GHZ state. Here, and throughout the manuscript, two 
states are termed "equivalent", if they differ only by a 
sequence of single qubit rotations, that is, local unitary 
operations. On the other hand, GHZ states belong to 
the class of graph states [19, 20] which can be gener- 
ated efficiently by applying two qubit operations (see 
Refs. [19, 38] and Chapter I). As a stabilizer quan- 
tum code the toric code ground states are always equiv- 
alent to a class of graph states. For the above example 
[Fig. 1 (A)] , the depicted equivalent graph state exhibits 
only connections (i.e. edges) between qubits which are 
vicinal within the toric code. However, this condition 
will break down for large setups: graph states, which 
are equivalent to the toric code ground states, will in 
general exhibit connections between qubits that are not 
nearest-neighbours, and therefore their experimental re- 



alization requires non-local interactions. Hence a funda- 
mental property of the toric code is the quasilocal nature 
of its stabilizers - which distinguishes it from the equiv- 
alent non-local graph states. Many applications of the 
toric code as topologically protected quantum memory 
rely on this local character. 

A minimal setup for topological matter should inherit 
this crucial property. In this manuscript, we present min- 
imal setups for the toric code, such that the equivalent 
graph states are non-local, to wit, they require connec- 
tions between qubits that are not nearest-neighbours. A 
major obstacle is that the equivalence between the toric 
code ground states and graph states is not unique, since 
(in general) there are multiple equivalent graph states. 
As a consequence, we have to exclude the existence of 
any equivalent graph state with only local connections. 
We therefore first derive a simple, geometric procedure to 
construct equivalent graph states from a given toric code 
setup. In a second step, all equivalent graph states for 
promising candidates are generated and tested on their 
property of non-local connections. 

The results for the most important setups are shown 
in Fig. 2: We find that for the common square lattice 
the minimal number of qubits is N = 16, see Fig. 2 (A), 
while for a triangular setup the minimal system reduces 
to N = 9 qubits, see Fig. 2 (B). 

Finally, we have to point out that our main interest lies 
in equivalent, that is local unitary (LU) equivalent graph 
states. However, computing LU-classes of graph states 
turns out to be exceedingly difficult whereas equivalence 
classes with respect to local Clifford (LC) operations can 
be computed in principle for any number of qubits algo- 
rithmically [20, 39-41], though not efficiently. For graph 
states it is known [19] that up to 7 qubits LU- and LC- 
classes coincide. Unfortunately, the conjecture that this 
equivalence holds for arbitrary numbers of qubits (see 
[42] and references therein) was disproved by means of a 
27-qubit counterexample [43]. However, it is believed [43] 
(but so far not proved) that this counterexample is the 
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smallest one. Since the systems considered in this letter 
are well below this threshold, we will deal with LC- in- 
stead of LU-equivalent graph states in order to simplify 
computations. 

This paper is organized as follows. In the first chapter 
(I) we summarize the mathematical and physical prereq- 
uisites used throughout this work: First, we outline the 
basics of the stabilizer formalism and give a short re- 
mark on different notions of "equivalent" states. Subse- 
quently, two important classes of stabilizer states, graph 
states and surface codes, are reviewed in more detail. In 
the second chapter (II) we present our first result, prov- 
ing that there is a simple geometrical procedure for the 
transformation of surface codes into LC-equivalent graph 
states. This procedure is illustraded by means of two ex- 
amples. In the last chapter (III) we present our second 
result: the smallest toric code setups for which no local 
graph state exists. Readers familiar with the basics of 
graph theory, graph states and surface codes may skip 
the preliminaries and proceed with Chapter II. Experi- 
mentalists primarily interested in the proposed non-local 
setups may omit the second chapter, as it is only of tech- 
nical interest for the setups in question. 



I. PRELIMINARIES 

The settings considered in this paper are composed of 
A qubits (or spins) and characterized by their coupling 
interactions. Thus the states are described by elements 
of the 2 Ar -dimcnsional Hilbert space H N := ®jL 1 Cf. 
The Pauli matrices acting on the Hilbert space of the ith 
qubit are denoted by of, a\, erf and, for the sake of con- 
sistency, cr? = 1 . As a basis of the 1-qubit Hilbert space 
rl 1 one chooses usually the computational basis {|0) , |1)} 
defined as eigenbasis of er z . Henceforth we shall use the 
convention a z |0) = |0) and a z |1) = — |I). 

In this paper we consider two distinguished classes of 
A-qubit states known as surface codes and graph states. 
Both are examples for the much larger and intensively 
studied class of stabilizer states. Thus we will start with 
a short survey of the stabilizer formalism. 

This formalism is widely used in quantum informa- 
tion theory and describes a certain class of N-qubit 
states, called stabilizer states, and their transformations 
by means of stabilizer groups. 

To this end define the Pauli group G\ := 
span {o~ x , cr^ , o~ z } operating on a I-qubit Hilbert space 
and generated by the Pauli matrices. The generaliza- 
tion for A qubits is defined as the A-fold tensor prod- 
uct Gn '■= of the 1-qubit Pauli group. Now 
let Q = {<7i} i=1 C Gn \ {—1} be a given set of in- 
dependent, pairwise commuting Pauli operators. Then 
S := span (J is called stabilizer of the protected subspace 
VS = {|$) G rl N |5|$) = |$)}. G is the generating set 
of iS and d = \Q\ =: rank S is called rank of the stabilizer. 
The dimension of the proteced space is dimVS ~ 2 N ~ d 



(for a proof and further reading see [44]), thus there are 
2 independent stabilizer states in VS. 

It proves advantageous to describe A-qubit states 
1 4') G VS by means of their stabilizing operators gen- 
erated by Q, for one has to track at most A genera- 
tors whereas the description of general A-qubit states 
requires an overall number of 2 N complex coordinates. 
In order to describe unitary transformations of stabi- 
lizer states by transformations of their stabilizer, one may 
use only unitaries that leave Gn invariant under conju- 
gation. In this context, one defines the Clifford group 
Ci := N u(2) (Gi) = {U e U (2) | UdW = Gi} as nor- 
malizer of the Pauli group in the 1-qubit unitary group 
U(2), and canonically extends this definition to the lo- 
cal Clifford group 

C l N : = <8>?=i C i ] operating on an A- 
qubit Hilbert space. If the transformations of a given 
stabilizer state are restricted to C l N (or more generally 
Cn-, see [45]), it is sufficient to transform the genera- 
tors Q — > CQC^ to describe the transformation G € C l N 
in the stabilizer formalism since CQC^ once more meets 
the requirements of a generating set. 

Whereas in the stabilizer framework local Clifford op- 
erations play a crucial role, one generally deals with a 
larger class of local operations, namely local unitaries 
U € U l N := f/(2)W. Since two arbitrary multi- 

qubit states cannot be transformed into one another by 
local Clifford operations or even local unitaries in general, 
the class of A-qubit states decomposes into equivalence 
classes with respect to restricted sets of transformations. 

To this end let and \*&') be two arbitrary A-qubit 
states. They are called equivalent under local unitary 
operations [LU-equivalent or equivalent for short) if there 
exists U £ U l N , such that |\fr) = U\H>')\ write |*) ~ LU 
I*!*'). Accordingly one calls them equivalent under local 
Clifford operations (LC-equivalent) if there exists G G 
C l N , such that |*) = G|*'); wri t c I*) ~lc I*')- 

Since C l N is a proper subgroup of U l N , it follows LC => 
LU but as C l N ^ U l N the converse cannot be true in gen- 
eral. However, due to the lack of counterexamples it was 
conjectured [42] that for the smaller class of stabilizer 
states LU- and LC-equivalence coincide: LC -O- LU. Al- 
though it was shown that this holds for a large class of 
stabilizer states [46, 47] the LU-LC-conjecture was dis- 
proved by Ji et al. in 2009 [43]. Nevertheless it was 
shown in Rcf. [19] that LU LC holds for graph states 
with up to 7 qubits. Furthermore it is believed [43] (but 
so far not proved) that the counterexample used in Rcf. 
[43] is the smallest one (A = 27). Recently the class 
of stabilizer states for which the LU-LC-conjecture holds 
was enlarged [48], including a large class of surface codes. 

After this brief survey of basic concepts of the stabi- 
lizer formalism and the notion of LC-equivalence we may 
now define a special class of stabilizer states, called graph 
states. 

Graph states were introduced by Hein et al. [19] in 
2004 and denote a class of A-qubit states that can be 
described by simple graphs. In this context it is shown 
[19] that the aforementioned LC-equivalence is closely 
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related to the concept of local complementations known 
in graph theory. To state the concepts of graph states 
more precisely, basics of graph theory arc required. Thus 
we start with a short review of the mathematics that 
may be skipped by the reader familiar with the concepts 
of graph theory. 

A finite undirected graph is a pair G = (V, E) where V 
denotes a finite set of vertices and EC {{a, b} \a,b £ V} 
a set of edges connecting these vertices. If there are nei- 
ther loops (i.e. edges connecting a vertex to itself) nor 
multiple edges (i.e. vertices connected by more than just 
one edge) G is called a simple graph. In contrast, if there 
are multiple edges (i.e. E is a multiset) but no loops 
[49] G is called a multigraph. In the following, simple 
graphs shall be denoted by upper-case letters G,H, . . . 
and multigraphs by calligraphic letters C,M,. . . . 

For a given graph G = (V, E) the functions V(G) and 
E(G) denote the sets of vertices and edges, i.e. V(G) = V 
and E(G) = E. The order of a graph G = (V, E) is the 
number of its vertices and denoted by \G\ := \V\. One 
naturally defines the neighbourhood of a vertex v £ V as 
the subset N v := {w £ V | {w,v} £ E} C V of vertices 
directly connected to v. For the sake of brevity we shall 
denote the complete graph over a vertex set V by (V) := 
(V, {{v 7 w} \ v 7 w £ V; v ^ id}), that is, the unique simple 
graph with edges connecting any pair of vertices. 

Graph theory deals not only with properties of graphs 
but also with transformations of the latter. A particu- 
larly interesting transformation in the context of graph 
states is known as local complementation. The local com- 
plement t v (G) of a simple graph G = (V,E) at vertex 
v £ V is formally defined as t v (G) := G + (N v ), where 
+ denotes the binary addition of edges. Practically the 
procedure described by t v reads as follows. Focussing ex- 
clusively on the neighbours of v one deletes all edges in 
E connecting these vertices and adds new edges between 
formerly unconnected neighbours. This yields t v (G). 

In conclusion we note that two graphs G = (V,E) 
and H = (V, F) over the same vertex set V are called 
equivalent under local complementations (LC- equivalent) 
if and only if there is a sequence of local complementa- 
tions t = t Vi o t V2 o • • • o t vn such that G = t(H)\ write 
G ~lc H. In the context of graph states, overloading 
"LC-equivalcnt" proves convenient and is therefore in- 
tended; this will become clear below. 

We may now give the formal definition of graph states 
and comment on their properties. To this end let G be 
an arbitrary simple graph with edges E and vertices V. 
Consider a system of N = \G\ qubits - each identified 
with a vertex v £ V - and define the (independent) op- 
erators 

K V G := o% J] al £ G N . (1) 

w£N v 

Then Q[G] = {Kq | v £ V} is a generating set of the graph 
state stabilizer S[G] = spanC?[G]. If there is no danger 
of confusion, we shall omit the dependence on G, viz. 
S = S[G]. Since there are N independent generators, the 



protected space VS[G] is one-dimensional. The state sta- 
bilized by S[G], henceforth characterized by G, is called 
graph state and denoted by |G). 

In this paper we shall define graph states in this sta- 
bilizer representation. However, there is an alternative 
representation which proposes a simple procedure for the 
construction of graph states in terms of controlled z-gates 
U ab = ® l fc + Vl>~ ® a z b . Here = ±[l v ± «tJ] 

denote the projectors operating on qubit v. Then it can 
be shown [19] that 

\g)= n ^i+) v(g) ( 2 ) 

{a,b}£E(G) 

where |+) = l/y/2 (|0) + |1)) and |+) V(G) denotes the 
iV-fold tensor product of states indexed by the vertices 
V(G). According to Eq. (2) the empty graph on N ver- 
tices V is represented by the product-state \+) V ■ One 
may add an edge between two vertices a,b £ V by apply- 
ing a controlled z-gate U a b on the corresponding qubits. 
This yields a handy procedure for the construction of any 
graph state |G) on N qubits - edge by edge - provided 

an entropy free initial state |+) V ^ G ' ) . 

Apart from their implementation, graph states are par- 
ticularly interesting multi-qubit states since several phys- 
ical properties can be related to graph theoretic proper- 
ties of the underlying graph [19]. One of these prop- 
erties is the aforementioned local Clifford equivalence: 
Let |G) and \H) be two graph states with corresponding 
graphs G = (E, V) and H = (F, V). Then they are LC- 
equivalent if and only if their graphs are connected by a 
sequence of local complementations. Formally this reads 
|G) — lc \H) G —lc H. The notion of LC-equivalent 
graphs was known in graph theory for a while and it 
was shown that there is an efficient algorithm to check 
whether two graphs are LC-equivalent [39, 40]. Thus the 
aforementioned relation between the LC-equivalence of 
graph states and graphs provides an efficient algorithm 
to decide whether two graph states belong to the same 
LC-class. 

For this reason the question arises whether this insight 
into the nature of LC-equivalence may be extended to a 
larger class of multi-qubit states. Indeed, it was shown 
[41, 50] that graph states are universal standard forms of 
stabilizer states. That is to say, for any A^-qubit stabi- 
lizer state | "J) there exists an LC-equivalent graph state 

LC 

|G*). Since the transformation |*) — > |G*) is effi- 
ciently computable, the problem of LC-equivalence for 
stabilizer states is reduced to the class of graph states - 
where it can be solved by means of mere graph theory. 

Besides graph states the stabilizer formalism is widely 
used for the description of quantum codes. A prominent 
example we are particularly interested in is the toric code 
model (tcm) and its generalization, the surface codes. 

In 2002 Kitaev proposed the qubit model named toric 
code [2] showing the inherent capability of quantum error 
correction on a physical level and a topologically ordered 
phase at zero temperature. Whereas the tcm is typically 
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Figure 3. (Color online) The tcm with qubits ( O ) located on 
the edges of a square lattice embedded into a 2-torus. The a x 
and er z Pauli operators are denoted by© and®, respectively. 
The intersections of X- and Z-type loop operators, responsible 
for their commutation relations, are highlighted by ©. The 
upper [left] border is identified with the dashed lower [right] 
border. 



defined on a square lattice embedded into a 2-torus, one 
may generalize this idea to arbitrary (multi-) graphs on 
orientable compact manifolds, called surface codes [2, 48, 
51]. _ 

Aiming at the definition of surface codes, we need some 
notation regarding graph embeddings. To this end let £ 
be a multigraph and E a closed, orientable 2-manifold. 
A graph embedding X = Xc,n of £ into E can be in- 
formally understood as a "drawing" of £ on the surface 
E without intersecting edges. In the following "embed- 
ding" has always to be read as 2-cell embedding, that is, 
a drawing on the surface, such that every face is home- 
omorphic to an open disc in R 2 [52]. The genus g[E] of 
E is its number of "holes" or "handles". The genus g[C] 
of a graph is defined as the least genus of surfaces, such 
that an embedding exists. Each embedding Xc.s induces 
a set of faces denoted by ¥x(C) or simply P(£). Given 
a graph embedding 3ic,Y, the dual graph £* = (V*,E*) 
is defined so that V* := ¥%(£) and two dual vertices 



fied with the edges of £ (cf. graph states!). 
star and plaquette operators 



A. 



and 



Define the 



(3) 
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for each vertex s £ V(£) and face p £ ¥%(£). 
Here we identify a face or plaquette p £ ¥%(£) 
with its set of peripheral edges and a vertex v £ 
V(£) with its adjacent edges. Then S = S[X] :— 
span {A s , Bp \ s £ V(£), p £ ¥x{£)} is the stabilizer of 
the surface code corresponding to the chosen embedding 
Its degeneracy is dim "PS = 4 9 = 2 2g in the case 
of closed boundary conditions [51] as one can show eas- 
ily using Eulers formula for 2-cell embeddings [53]. A 
possible basis of the protected space PS is defined and 
manipulated by means of loop operators 



Z[C k 



of, X k = X[C%} 



(4) 



i£C k 



'ect 



where the simple (co-)cycles C k C £ and C £* are cho- 
sen homologously distinct for different k and non-trivial 
(see Fig. 3). In algebraic topology, two paths on a com- 
pact, orientable 2-manifold E are considered homologous 
if their union is the boundary of a surface embedded into 
E. In our case, two cycles on the embedded graph £ are 
homologous iff their union is the boundary of a union of 
faces. This defines an equivalence relation ~ on the group 
of all cycles £ (£). The quotient Hi := £ (£) is called 
first homology group and contains homology classes [C] of 
cycles C £ £(£). Restricted to the protected space VS, 
the effect of a loop operators Z [C] (X [C*]) depends solely 
on the homology class of C (C*) since they act trivially if 
the (co-)cycle is homologeously trivial, that is, C (C*) is 
the boundary of a union of faces (stars, i.e. dual faces). 
Therefore, by a slight abuse of notation, we may regard 
Z [C] cither as a mapping Z[ ■ ] with domain £ (£), or as 
an operator Z- with domain H\. 

If the indices k are chosen properly, the loop operators 
fulfil the algebra {Z k ,X k } = = [Z k ,Xi] for k ^ I. Then 
{|v) eVS\Z k |v) = v k jv) , 1 < k < 2g} is a basis of VS, 
where the components of v £ { — 1, l} 29 arc called topo- 
logical quantum numbers and represent the logical qubits 
encoded in the protected space. The topological quantum 
numbers arc flipped by the 2g distinct X k loop operators 
via X k \vi, . . . ,Vk, ■ ■ ■ ,v 2g ) = \vi, . . . , -v k , . ■ ■ ,v 2g ) due 
to the aforementioned (anti-)commutation relations. 

To conclude this section, we note that in the special 



£ V* are adjacent, i.e. {v^v^} £ E* , iff there case of a square lattice £ and E = T 2 a 2-torus, as de- 



is a common edge in E(£) to vl and u| as faces. A (co- 
)cycle C C E(£) (C* C E(£*)) on £ (£*) is a set of edges, 
such that each vertex is adjacent to an even number of 
these edges. Arc there at most two of these edges joining 
at each vertex and is C (C*) connected within the graph, 
then it is called a simple (co-) cycle, namely, a single cycle 
without self-intersections. 

For the definition of surface codes, consider a multi- 
graph £ 2-cell embedded via %c.s into a surface E of 
genus g[E] > g[C] and a system of |E(£)| qubits identi- 



pictcd in Fig. 3, this construction yields the tcm [2]. 



II. MAPPING SURFACE CODES TO GRAPH 
STATES 

A. Results 

In the following, let |v) be the state of a surface 
code with topological quantum numbers v, correspond- 
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Figure 4. The graph on the left hand side is simple and not 
a tree since it features cycles. Deleting edges that support a 
cycle (illustrated by dotted lines) yields two non-isomorphic 
spanning trees. 

ing stabilizer S[X] and embedded multigraph £. \G) de- 
notes a graph state with simple graph G and stabilizer 
S[G]. Since every stabilizer state is LC-equivalent to 
some graph state, there exists a local Clifford operation 
C € C l N such that \G) = C|v). We are now interested 
in finding C and S[G] for a given surface code S[X] in 
state |v). The binary algorithm developed by Van den 
Nest et al. [41] allows one to perform this computation 
algorithmically. Here we show that in the case of sur- 
face codes there is a purely geometrical, graph theoretic 
transformation rule. 

To this end a further notion of graph theory is required. 
Let £ and £' be (multi-)graphs. £' is called a spanning 
tree of £ if it is simple, connected, devoid of cycles and 
V [£] = V [£'] as well as E [£'] C E [£]. That is to say, 
spanning trees are distinguished subgraphs without cy- 
cles that cover all vertices of the original graph. However, 
they are not unique in general. Spanning trees of a given 
graph may be constructed by deleting a set £ C E[£] of 
edges such that there is no cycle left and the graph re- 
mains connected. For an example of this procedure see 
Fig. 4. Let us now introduce a mapping that converts 
multigraphs to a particular class of simple graphs which 
proves crucial for the LC-transformation of surface codes 
to graph states. 

Definition 1 (ipc> ■ £ n> G). Let £ be a multigraph 
and CJ one of its spanning trees. Then G = ipe' (£) * s a 
simple graph defined by the following construction: 

Set V(G) := E(£) as the vertex set of G. For r, s £ 
V(G) the set {r, s} is an edge of G if and only if there is 
a path C pq from p to q on £' (p,q £ V (£') = V (£) ) such 
that r = {p, q} and s £ C pq . 

The physical relevance of this mapping from multi- 
graphs to simple graphs is given by the following 

Theorem 1 (LC-equi valence). Given the surface code 
state |v) on £ and any spanning tree £'. Then the graph 
state \G) = |</?£'(£)) £ H N is LC-equivalent to |v). 

As a final remark note that the combination of the 
LC-rule for graph states [19] and Theorem 1 yields the 
purely graph theoretic statement that for arbitrary span- 
ning trees C[ and C' 2 it holds <pc' {C) <Pc (£)■ F° r 
the sake of consistency, we derived a direct graph theo- 
retic proof for this very statement (which shall be omitted 
here). 




Figure 5. (Color online) Two simple examples of surface 
codes and LC-equivalent graph states. Product states (A) 
are mapped to disjoint graphs and GHZ states (B) are trans- 
formed to star graphs. A description of the graphical notation 
is given in Subsection II B. 



B. Examples 



Here we discuss two examples, depicted in Fig. 5, in 
order to demonstrate the construction implied by Def. 
1 and show the consistency with well known facts. The 
graphical notation is as follows: Qubits are denoted by o 
whereas • denotes the vertices of the surface code, and 
thus the position of star operators A s . The edges of the 
surface code are denoted by thin, continuous and dashed 
lines connecting the vertices; the "deleted" edges £ are 
drawn with dashed, the spanning tree with continuous 
lines. Faces are highlighted grey and symbolize plaquette 
operators B p . Hadamard gates transforming qubits on 
"deleted" edges (see proof of Theorem 1 in subsection 
II D) are symbolized by®. LC-equivalent graph states 
are drawn with bold red edges connecting the qubits. 

Example 1 (Product states). It is easy to see that pro- 
tected states of two one-point connected surface code 
systems (Fig. 5 (A)) are product states, symbolically 
I CO) = |G>) ® I©). As a consequence of the one-point 
connection, any LC-equivalent graph state turns out to 
be disconnected - independent of the spanning tree used 
for the construction. Since disconnected graphs denote 
the tensor product of graph states [20] , we confirmed the 
statement |CO) = |G>) ® |©) in a purely geometric way, 
for local unitaries cannot alter the degree of entangle- 
ment. 

Example 2 (GHZ states). One can show eas- 
ily that the protected state of a single plaquette 
bounded by N qubits is LC-equivalent to the A^-qubit 
GHZ state |GHZtv) since it has the form |0) = 

l/y/2 (\+)® N + I-)®*) ■ Consider the hexagonal pla- 
quette depicted in Fig. 5 (B). Constructing spanning 
trees is easy since "deleting" an arbitrary edge yields one. 
Consequently any graph state constructed this way must 
be a star graph. Indeed, in [20] it is shown that star 
graphs (and complete graphs) are LC-equivalent to GHZ 
states. 
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C. Comments on the structure of tp£j(£) 

The geometrically constructed graph ipc (£) exhibits 
an interesting graph theoretical structure, which in turn 
furnishes insight into the entanglement structure of the 
corresponding surface code state. To this end, we intro- 
duce the cycle space £ (£) and its dual, the cut space 
as linear subspaces of the edge space £(£); the 
latter being the power set of edges E (£) equipped with 
binary addition and multiplication. The cycle space C (£) 
is defined as set of all cycles in £ (as sets of edges) 
equipped with the binary structure inherited by £(£). 
A cut of £ is a bipartition of its vertices; the cut-set of 
a given cut is the set of edges in £ connecting these two 
disjoint sets of vertices (often the cut-set is just called 
cut). The cut space D (£) is then defined as set of all 
cut-sets on £ equipped with the usual binary structure. 
Bases of both subspaces may be constructed by means of 
an arbitrary spanning tree £' as follows. For each edge 
e = {Pi q} S E (£) \ E (£') choose the cycle C e := C pq + e 
where C pq denotes the unique path on the spanning tree 
connecting p and q. C e is called fundamental cycle on £ 
with respect to £' and one can show easily [54] that the 
set of all fundamental cycles {C e | e £ E (£) \ E (£')} is a 
basis of the cycle space. Conversely, each edge e' G E (£') 
induces a bipartition VJ: U¥j = V (£) of the vertices, 
according to the component of the spanning tree they 
belong to if e' is deleted. The corresponding cut-set T> e i 
is called fundamental cut of e' with respect to £'. Anal- 
ogously the set of all fundamental cuts {V e > \e' G E (£')} 
yields a basis of the cut space. 

Within this graph theoretical framework, the structure 
of ipc' (£) reads as follows. First, we notice that (fc (£) is 
bipartite, that is, there is a bipartition of its vertices, so 
that edges interconnect vertices of different classes only. 
This bipartition is given by the set of vertices which be- 
longed to the spanning tree £' as edges of £, and its corn- 




Figure 6. (Color online) Interpretation of the bipartite graph 
G = <fic'{£) as given in Def. 1. Each vertex e identified with 
a deleted edge in £ (#) encodes a fundamental cycle C e on £ 
with respect to £' via its neighbours in G, viz. C e = N e U 
{e}. Conversely, each vertex e' identified with a spanning tree 
edge (■) encodes a fundamental cut T) e i on C with respect to 
C! via its neighbours in G, viz. T> e > = N e / U {e'}. This is 
a consequence of the duality between fundamental cuts and 
fundamental cycles. 



plement. Therefore wc may draw tf£i (£) schematically 
as depicted in Fig. 6. It is now straightforward to see that 
the neighbourhood N e (N e > ) of vertices e G E (£) \ E (£') 
(e' G E (£')) on tpc(C) joined with e (e') equals the fun- 
damental cycle C e (cut T> e >) on £ with respect to £'. 
We conclude that ipc(C) encodes the whole cycle space 
£ (£), and the cut space T) (£) as well. In Fig. 6 the high- 
lighted subgraphs represent a fundamental cycle C e and 
a fundamental cut V e i indexed by e and e', respectively. 
That both, <£(£) and D (£), are encoded in the same 
graph is not a conincidence but rather a necessity since 
it is well known that there is a duality between funda- 
mental cycles and fundamental cuts [54] ; ipc' (£) is a ge- 
ometrical illustration of this duality. Since entanglement 
within the framework of graph states is characterised by 
edges between the qubits/ vertices, the LC-equivalence of 
[</?£'(£)) to a given surface code ground state on £ illus- 
trates that the entanglement pattern of a surface code is 
basically determined by its cycle and cut space structure. 

D. Proof of Theorem 1 

In this subsection the proof of Theorem 1 is presented 
in detail. Readers primarily interested in the minimal 
non-local toric code settings may skip this rather techni- 
cal part and proceed with Section III. 

Proof of Theorem 1. First note that the basis states |v) 
of the protected space are LC-equivalent since the topo- 
logical quantum numbers Vi are controlled by loop opera- 
tors Xk G C N . Therefore, from the viewpoint of LU/LC- 
classification, w.l.o.g. assume v = (1,...,1). It follows 
that for each cycle in £ there is a z-type loop operator 
in the stabilizer S[X] (Fig. 7 (A)). Let £' be an arbitrary 
spanning tree and £ the set of edges such that £' = £ — £ 
(Fig. 7 (B)). Apply Hadamard transformations H e G G N 
to the qubits identified with edges e G £. Formally this 
transformation reads |v) i-> H |v) where H £ = ® e ^sH e . 
We now argue that H £ |v) = |</?£<(£)) by showing that 

the new stabilizer S' := H £ S[X] (H £ ) is the graph state 
stabilizer described by ip £<(£). 

In the stabilizer picture, a Hadamard transformation is 
represented by the substitution a x <-> cr z . Consider one 
of these transformed qubits (® in Fig. 7 (C)) and its edge 
e = {p, <l} € E(£) (p, q G V(£)). The adjacent vertices p 
and q belong to both £ and £'. Since £' is a tree, there 
exists one and only one path C pq on £' connecting p and 

q ( in Fig. 7 (C)). Let u z e ® Z [Cpq] be a combination 

of plaquctte operators in S[X] (note that e-\-C pq is a cycle 
in £). This stabilizer operator is transformed under the 
Hadamard transformation as follows 

K { * ] = H e a z e <g> Z [C pq ] Hi = a x e ® (g) of . (5) 

Obviously this is a graph state generator of <S [</?£/(£)] 
since in <pc' (£) vertex e is connected to all edges in £ 
that belong to C pq (Fig. 7 (D)). This argument holds for 
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Figure 7. (Color online) Construction of G = ip^'iA) as given in Def. 1 and used in the proof of Theorem 1. Given a surface 
code (A) with arbitrary spanning tree (B). Consider one of the (deleted) edges e = {p, q} on which a Hadamard rotation was 
performed (denoted by®) and its unique path C pq (denoted by - - - ) on the spanning tree (C). The former z-type loop operator 
along Cpq + e acts as a x on the qubit © (due to the Hadamard rotation) and as a z on the qubits denoted by ®. This is a graph 
state stabilizer of the graph drawn with (D). On the other hand, consider one of the eges e' on the spanning tree (E). The 
induced bipartition N% U N% of the vertices is denoted by A and ■. The product of all (former) star operators labeled by one 
of these sets acts nontrivially as a x on the qubit © and as a z on the qubits denoted by ® (due to the Hadamard rotations). 
This is a graph state stabilizer of the graph drawn with in (E). 



any edge in £. Therefore we find that the transformed 
stabilizer S 1 contains at least the graph state generators 
of ip£/(£) for vertices in £ . It remains to show that the 
graph state generators for vertices e' £ E(C') can be 
found in S' as well. 

To this end consider an arbitrary qubit e' G E (£') in 
the spanning tree (© in Fig. 7 (E)). The corresponding 
edge of the spanning tree cuts the latter into two parts. 
This construction yields a bipartition of the vertex set 
V (£') = Vf,' U Vjf (red A and blue ■ in Fig. 7 (E)). Thus 
there are two types of edges in £ : The first type connects 
vertices of the same partition and the second type inter- 
connects both sets and . By definition of ipc (C) 
a vertex in £ is connected to e' iff it is a type-two vertex. 

We claim that the following operator in S' is the graph 
state generator for e' 

K g^ = II A 's whcrc A' 3 =H £ A a (H s y . (6) 

To see this, note that risGV e ' may affect a qubit / in 
three different ways: 

First, if / belongs to an edge with adjacent vertices of 
the same partition, nothing happens since o~ z a z = 1 = 
a x a x . The second case occurs if the adjacent vertices 
belong to different subsets and f ^ e' (® in Fig. 7 (E)). 
Since this case can only apply if / G £, a a z opera- 
tion is performed on this qubit. The last case resembles 
the previous one, except for the fact that / = e! . Thus 

riseV e ' acts as ° x on /■ 

To sum up: we found that Kq acts as o~ x on e , 
as a z on all vertices / such that {e',/} <G E(ipc(£)), 
and trivially otherwise. Therefore we showed that 
S [<p£'(£)] < S' . To show S [<pc'(C)] — S' we note that 
rankS [<pc'(C)] = N = rankS [X] = rankS'. 



III. MINIMAL NON-LOCAL TORIC CODE 
INSTANCES 

Setting the scene. In this section, we present mini- 
mal instances for the toric code on square and triangular 
lattices which are not LC-equivalent to any local graph 
state. We already provided a brief explanation and moti- 
vation for the notion of "local" in the introductory para- 
graphs. At this point, we shall give a more formal def- 
inition and relate the idea of locality to the toric code 
hamiltonian. 

The latter is usually written as 

H tcm = -J A A S -J B b p (7) 

s£V(£) pGP(£) 

where J a and Jb define the excitation energies for elec- 
tric and magnetic charges, respectively. Obviously the 
aforementioned protected space PS is also the ground 
state space of the above hamiltonian and the orthogonal 
states |v) constitute a basis of the latter. In the prelim- 
inaries the toric code was described as a quantum code 
within the stabilizer framework. However, considering 
the hamiltonian (7), the toric code gains significance as 
an actual physical system that may be implemented di- 
rectly or by means of a quantum simulator. Crucial for 
the physical relevance of H tC m is its quasi-locality, that 
is, only quasi-local interactions - represented by star and 
plaquette operators - are required. Therefore, from a 
physical point of view, the star and plaquette operators 
are distiguished by their quasi-locality, whereas from a 
mathematical point of view all generating sets of S are 
equivalent. 

Consequently, the set of star and plaquette operators 
gives rise to a physically motiviatcd adjacency relation 
on the set of qubits - namely, two qubits p and q arc 
vicinal if both belong to the support of at least one star 
or plaquette operator, write p h c q. On the other hand, 
the relevance of graph states in experimental physics is 
due to their scalability which in turn relies on their suc- 
cessive implementation by means of 2-qubit operations 
(see Eq. (2)). Therefore it seems appropriate to call two 
qubits vicinal in a graph state if they are connected by 
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an edge, write pl gs q. The introduced symmetric relations 
I encode physically motivated adjacencies with respect to 
a particular class of multipartite entangled states. 

After introducing the notion of vicinity, we may now 
ask the following question: Given the ground state |v) of 
a toric code, is there an LU-equivalent graph state |Gi oc ), 
such that for all pairs of qubits pl gs q p he 1 holds. 
Provided such a graph state exists, we may apply those 
implementation techniques for graph states which gave 
rise to l gs (within the setup defined by the toric code) in 
order to generate its ground states. Conversely, if there 
is no such state, the considered toric code setup yields 
ground states which are fundamentally different from its 
equivalent graph states, that is, the toric code ground 
states cannot be implemented by means of graph states 
without breaking the spatial structure of the toric code 
setup. At this point, the quasi-local property of -fftcm 
becomes relevant. Thus we call toric code setups that ex- 
hibit at least one graph state |Gi oc ) local and G\ oc (|Gi oc )) 
itself a local graph (state) with respect to the considered 
toric code setup. Hence a toric code is local if one can 
draw an LU-equivalent graph state without "long-range" 
edges, that is, connecting qubits which cannot interact 
via star or plaquette operators in the toric code regime. 
E.g. see Fig. 2 for non-local and Fig. 5 for local instances 
of toric code systems. 

Small toric code systems tend to be local, e.g. the 
single toric code plaquette in Fig. 1 (A) is represented 
by a local star graph. However, this property is likely to 
break down for more complex setups. This motivates our 
search for the smallest non-local toric code systems. In 
doing so, we restrained the analysis to small systems on 
square and triangular lattices, since experimental imple- 
mentations are likely to use such tesselations. 



A. Methods 

Our approach is as follows. Due to the restriction to 
a particular tesselation (let us consider the square lat- 
tice exemplarily) we are looking for connected subsets 
of squares on the square lattice. To this end we may 
consider only setups with edge-connected squares since 
vertex-connected structures yield separable ground states 
or (if there is a loop through the vertices) degenerate 
ground states; the latter requiring a reduction to a non- 
degenerate setup by "filling" the loop, anyway [55]. Such 
structures are known in combinatorics as polyominos (on 
the square lattice) and polyiamonds (on the triangular 
lattice). The polyominos for up to n = 5 squares are 
depicted in Fig. 8. Application of Theorem 1 yields lo- 
cal graph states for all depicted setups devoid of "inner" 
squares - provided a suitable spanning tree is used for the 
construction. Actually, by means of Theorem 1 this can 
be worked out mentally. Subsequently there is a smallest 
setup left, such that the geometric transformation yields 
only non-local graphs. For the square lattice this would 
be the +-shaped pentomino (n = 5) with dark coloured 



□ 



Figure 8. (Color online) The free polyominos for n = 
1, 2, 3, 4, 5. Faces with at least one peripheral edge are light- 
coloured. The first polyomino featuring a periphery-free cell 
(dark blue) is the +-shaped pentomino. It is the smallest 
configuration on the square lattice, such that the geometical 
construction yields only non-local graph states. 



inner square, as depicted in Fig. 8. 

However, the class of geometrically constructed graph 
states is generally a proper subset of the whole LC-class 
containing all LC-equivalent graph states [56]. To fix 
this issue, we computed the whole LC-class and searched 
for local representatives. To this end, an algorithm de- 
scribed in [20] and originally derived in [39, 40] was em- 
ployed, according to which two graphs G and G' with 
adjacency matrices [57] T and T' are equivalent under 
local complementations iff there exist diagonal matrices 
A,B,C,D e xN , such that the non-linear condition 
AD + BC = E and the system of linear equations 

(rB + D) r' + (rA + c) = o (8) 

are satisfied over the Galois field F2 (E denotes the 
identity matrix). This condition was implemented using 
Mathematica and facilitated the generation of the whole 
LC class [G]lc f° r a given graph G. However, this al- 
gorithm is inefficient in the sense that computation time 
increases exponentially with the number of vertices. To 
our knowledge, no efficient (i.e. polynomial) algorithm 
has yet been found. Unfortunately, the exponential scal- 
ing behaviour renders the LC-class generation for the +- 
shaped setup (i.e. N — 16 vertices/qubits) unfeasible 
on a regular personal computer, whereas the triangu- 
lar tetriamond depicted in Fig. 2 (B) (i.e. N = 9 ver- 
tices/qubits) requires only few minutes to generate the 
whole LC-class and seek for local graphs. To deal with 
this (rather technical) problem, we conceived a method 
to deduce the non-locality of certain toric code setups 
from smaller ones. The details concerning this method 
are technical and purely graph theoretic. We therefore 
omit them at this point and refer the interested reader 
to Appendix A. 
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B. Results 



As described above, the smallest setup on the square 
lattice without geometrically producable local graph 
state is the +-shaped pentomino (N = 16), as depicted in 
Fig. 2 (A). The aforementioned reduction method yields 
the smaller (N — 8) setup depicted in Fig. 2 (C). We 
proved the non-locality of the latter by generating the 
whole LC-class and probing for local graphs. Since the 
result was negative, we conclude that the system depicted 
in Fig. 2 (C) as well as the 5-plaquette setup in Fig. 2 (A) 
are non-local. Therefore, on the square lattice, the +- 
shaped pentomino with N — 16 qubits is the smallest 
setup without (LC-)equivalent local graph state. 



On the triangular lattice the triangle shaped tetria- 
mond in Fig. 2 (B) is the smallest one yielding only non- 
local, geometrically constructed graph states. Further- 
more, the generation of all LC-equivalent graphs yields 
only non-local ones. Thus, on the triangular lattice, the 
triangle-shaped tetriamond with N = 9 qubits is the 
smallest setup without (LC-)equivalent local graph state. 
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Appendix A: Reduction of surface code systems 

1. Structure of e-symmetric LC-classes 

In this section, some properties of LC-classes of simple 
graphs are presented that may be used to decide whether cer- 
tain surface code systems are local. The following statements 
are purely graph theoretic. 

Lemma 1 (e-Symmetry). Let [G]lc be the LC-class of a sim- 
ple graph G such that there exist representatives with at least 
one leaf. If a and b denote the outer and inner vertex of a 
leaf, then the whole LC-class is invariant under the exchange 
of a and b. We denote this symmetry by 

eab ■ <£{V) -> £(V) (Al) 

and call [G]lc e-symmetric with respect to a and b. Formally 
this reads e a b[G]i,c = [G]lc- 

Here £(V) denotes the edge space over the vertex set V, 
that is, the set of all simple graphs with vertices V equipped 
with the binary addition of edges. 



Proof. Let G a b be a representative with a leaf and a be the 
outer vertex. Then it is easy to verify that the two local 
complement at ions 

Gba = TaTbGab 

yield a graph Gba such that b is the outer vertex and all edges 
that connected b in G a b with vertices V \ {a, 6} now emanate 
from a. In other words: The graph Gba is obtained from 
G a b by exchanging the vertices a and b. Hence the symmetry 
£abG a b = Gba is realized by two local complementations. Now 
consider an arbitrary element H a b G [G]lc- Hab is reachable 
from Gab by subsequent application of local complementa- 
tions. If r(a,b) denotes this chain of LC-operations, then let 
r(b,a) be the chain obtained by the substitution r a -H- Tf,. 
Then H ba = r(b,a)G b a = r(b,a)e a bG a b = T(b,a)T a T b G a b dif- 
fers from Hab by a permutation of a and b; we write Hba = 
EabHab G [G]lc- Thus we find e„ 6 [G] L c C e a b[G]hc C [G]lc 
and since e' 2 ah = 1 it follows e a b[G]-Lc = [G]lc- ■ 

The symmetry originally derived for the special leaf-graphs 
Gab and Gba "spreads" through the whole equivalence class 
and leads to some useful properties regarding LC-classes rep- 
resented by such graphs. Before analysing this structure, let 
us prove the following statement. 

Lemma 2 (Deletion of leafs). Let G ab £ £(V) be a simple 
graph with at least one leaf; a denotes the outer and b the 
inner vertex. Consider the subgroup 

Ta := span {n : £(V) -)• &{V) \ieV\ {a}} (A2) 

of all local complementations on vertices V \ {a}. Then it 
holds 

Vr eT a : rGab - a = t (Gab - a) . (A3) 

In words: Deleting the leaf (a and its adjacent edge) from 
Gab and transforming this graph by a sequence of local com- 
plementations t acting on vertices V \ {a} yields the same 
graph as if these operations were performed on G a b (without 
deleting the leaf) and deleting a with all its adjacent edges 
(which may be more than one) subsequently. 

Proof. It is easy to see that all changes regarding edges be- 
tween vertices in V \{a} occur in the same way, whether there 
is an edge to a or not. Therefore the only differences between 
rG ab and r (G a b — a) affect edges from vertices in V \ {a} to 
a itself. Since these edges are deleted (rGab — a), there is no 
difference left and it follows rG ab — a = t (G a b — a). ■ 

Lemma 2 allows statements about the LC-class [G a b — o]lc- 
More precisely: All graphs in [G a b]LC that are reachable from 
Gab by concatenations of local complementations in T a yield 
graphs in [G ab — o\lc by deleting a and its adjacent edges. If 
one can show that all graphs in [G a f,]LC may be linked to each 
other by operations in such subgroups, this yields a procedure 
to compute LC-classes of special subgraphs. 

This idea leads to the following proposition for simple 
graphs with leafs: 

Proposition 1 (e-symmetric LC-classes). Let G ab be a sim- 
ple graph with at least one leaf where a and b denote the outer 
and inner vertex, respectively. 

Then for all H 6 [G a f,]LC there exists H' a G [G a b — h]lc 
or H' b e [G ba - &]lc such that H - a = H' a or H - b = H' b , 
respectively. 
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In other words: If an LC-class features a leaf-graph repre- 
sentative with outer and inner vertex a and b, then deleting of 
one of the latter (in some cases it does not matter which one 
is deleted) yields a graph in [G a b — o]lc or [Gba — 6]lc- Con- 
sequently, any graph in [G &]lc = [Gb a ]hC features at least 
one subgraph in [G a b — «]lc or [Gba — 6]lc- 

Proof. Let G a b be a representative with outer and inner vertex 
a and b. Then one can partition the LC-class [G a b]hC in the 
following way: 

• All graphs with a leaf belong to A iff a is the outer and 
b the inner vertex. Therefore G a b £ A. 

• All graphs with a leaf belong to B iff b is the outer and 
a the inner vertex. Therefore Gba £ B. 

• All graphs H where {a, b} G E(H) and for all v G V(H)\ 
{a, b} it holds 

{a,v} G K(H) {b,v} G E(H) 

belong to C. Due to the second statement, these graphs 
are totally symmetric with respect to a and 6, e a bH = 
H, that is. 

• All graphs H where {a, b} $ E(H) and for all v G V(H)\ 
{a, 6} it holds 

{a,v} G K(H) <=> {b,v} G E(H) 

belong to D. These graphs are totally symmetric; the 
only difference to type-C graphs is the missing edge 
between a and b. 




Figure 9. (Color online) Diagrammatic sketch showing the 
structure of e-symmetric LC-classes as used in the proof of 
Proposition 1. Arrows denote different classes of transfor- 
mations according to their color: 7-operations (magenta), 
^-operations (green), a-operations (blue) and fa-operations 
(red). A detailed description is given in the text. 

These four types of graphs are depicted in Fig. 9. Note 
that the number of edges between a or b and the vertices 
V(G) \ {a, b} (denoted by light grey discs) is not constant 
within one class in general. That the described classification 
of graphs is a partition, i.e. [GoJlc = AO B Cl CO D has to 
be shown. To this end, we introduce a classification of local 
complementations t v G T which depends on the graph they 
operate on. First, we call r a and Tb an a- and b-operation, 
respectively. Secondly, local complementations r v (acting on 
G) where v G V(G) \ {a, 6} and {a,v},{b,v} £ E(G) are 



called -^-operations and LC-operations t v (acting on G) where 
v G V(G) \ {a, b} and {a, v} G E(G) or {6, v} G E(G) are 
called ^-operations. Obviously every local complementation 
t v , v G V(G), belongs (depending on the graph it is applied 
to) to one of these four classes. Let us examine how these 
four classes of local complementations and the four classes of 
graphs introduced above collude: 

• First, note that 7-operations cannot alter the class a 
graph belongs to (A,B,C or D) since all edges used 
for the classification are unaffected by 7-operations. In 
Fig. 9 this is denoted by magenta arrows pointing from 
each set into itself. 

• Clearly, a-operations act trivially on A. Applied to 
an element of G, it deletes all edges between b and 
V(G) \ {a, b} and yields a graph that belongs to B. 
Since rf = 1, this is true for the opposite direction 
as well. Finally, graphs that belong to D cannot leave 
this set since no edge between a and b is created. These 
connections are denoted by blue arrows in Fig. 9. 

• Due to the symmetry with respect to a and b, all state- 
ments for T a hold for tj, (i.e. b-operations) as well. 
These relations are denoted by red arrows in Fig. 9. 
Note that r a equals t d if they are applied to elements 
in D. 

• ^-operations act on vertices denoted by dark grey cir- 
cles in Fig. 9 (vertices that are connected to a or 6 or 
both of them). Obviously such operations cannot ma- 
nipulate the leaf. Therefore ^-operations cannot leave 
the sets A and B (but they may alter the number of 
edges connecting a or b to vertices in V(G) \ {a, b}). 
Acting on elements in G, a ^-operation deletes the edge 
between a and b, hence converting type-G graphs to 
type-D graphs and vice versa. These relations are de- 
noted by green arrows in Fig. 9. 

Consider the following situation: Starting from G a b G A 
one tries to generate as much elements in [G a i>]LC as possi- 
ble by using local complementations in Ta- Therefore one is 
allowed to use b- 7- and ^-operations to generate graphs. In 
terms of Fig. 9 the "path" cannot use the blue arrows. Note 
that the only a-operation that cannot be omitted connects G 
and B (in A it acts trivially and in D it may be replaced by a 
fe-operation). We call graphs in this orbit T a G a b a-accessible. 
In conclusion it holds T a G a b Q [G]lc\-B. Analogue arguments 
yield TbGba C [G]lc \ A for b-accessible graphs starting from 
Gba £ B. 

We are now going to show that the more stricter relations 
TaGab = [G] LC \ B and TbGba = [G]lc \ A hold. To this 
end, consider an arbitrary "path" jumping from class to class, 
w.l.o.g. starting at G a b and ending at H G [G]lc\B, formally 
described by 

H = Y[ TiGab where n G T 

i 

as depicted in Fig. 10 (A). In general there will be some a- 
operations in this chain. Since neither G a b nor H are in B 
there is an even number of them (in Fig. 10 (A) there are 
two, denoted by blue jumps). We are restricted to non-a- 
operations in Ta, though. We can adhere to this restriction 
by "swapping" the a-jumps and the corresponding 7- and £- 
operations in B to &-jumps and the same 7- and ^-operations 
in A. This is shown in Fig. 10 (B). 
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Figure 10. (Color online) An arbitrary path from G a b (red 
dot) to H (green dot) is depicted in (A). The "forbidden" 
(blue) paths using a-operations are swapped to A (red jumps) 
using 6-operations instead, see (B). Now the whole path from 
Gab to H is a-accessible. A detailed description is given in 
the text. 



It is easy to see that jumping to A by 6-operations in- 
stead and subsequently applying the same 7- and ^-operations 
yields a graph in A that resembles the graph obtained in B 
before jumping back to C with the only difference that a and 
b are permuted. However, since the second application of t d 
yields a totally symmetric graph in C, this graph has to be 
the same one as formerly obtained by the second r a operation. 

Consequently we found a procedure to reach any graph 
H G [G]lc \ B by local complementations in T a from G a b 
and, due to symmetry, each H' G [G]lc \ A by exclusive 
application of Tb operations from Gba [58]. Thus we showed 
%Gba = [G] LC \ A and T a G ab = [G] hc \ B. 

To conclude the proof, recall Lemma 2. Consider H G 
[G]lc- Assume without loss of generality that H ^ B (the 
following argument holds for H ^ A analogously). Then we 
showed that H G T a G a b where G a b is a leaf-graph with outer 
vertex a. Choose r £ T a such that H = rG a b- According to 
Lemma 2 it follows that H — a = rG a b ~ a = r (G a b — o). Set 
H' a :— t (Gab — a) G [G a b — aJlc and it follows H — a = H' a G 
[Gab — «]lc. Which was to be proven. 



2. Application to the locality problem 

In order to apply the graph theoretic knowledge about e- 
symmetric LC-classes introduced above, we have to define in 
graph theoretic terms the notion of "local surface codes" and 
"local graph states" . 



Definition 2 (Adjacency relation, Local graph). Let V be 
an arbitrary finite set. An adjacency relation A on V 
is an irreflexive and symmetric relation, i.e. there is no 
v G V such that v ~a v and for each pair v,w G V it holds 
v ~a w =r* w ~a v. An adjacency relation is represented 
by an adjacency matrix A where A V:W = 1 if v ~a vj and 
A VjW = otherwise. Therefore A defines a simple graph Ga 
with adjacency matrix A. A simple graph Gi oc is called local 
with respect to A if it is a subgraph o/Ga. 

Given a surface code £ (more precisely: a 2-cell embedding 
Xc.s of £ in E), the stabilizer generators A s and B v induce 
an adjacency relation A on the set of edges (or qubits) E(£) as 
follows. Two qubits e, / G E(£) are adjacent, e ~a /, if and 
only if there is a common face p G P(£) or a common vertex 
s G V(£) such that e, / G p or e, / G s, respectively. Note 
that by this definition of A, local graphs G\ oc in the sense of 
Def. 2 are the very graphs describing local graph states |Gi oc ) 
according to the definition given in Section III. 

Since we want to reduce surface code setups to smaller ones 
by removing edges/qubits, a method to compare two setups 
with respect to their adjacency relations is needed. Therefore 
we introduce the notion of strictness. 

Definition 3 (Strictness). Let £1 and £2 be two surface codes 
with edges (qubits) £(£2) C E(£i). Furthermore denote by 
Ai and A2 the adjacency relations as defined above. They 
define simple graphs Gh t and Ga 2 with vertices E(£i) and 
£(£2), respectively. Ai is called stricter than A2 (write Ai )p 
A 2 ) if Gax [E(£ 2 )] C Ga 2 . 

Here Gai [E(£2)] denotes the subgraph of Ga ± on the ver- 
tices £(£2) C E(£i). Definition 3 and Proposition 1 enable 
us to state and prove the following theorem, which is the main 
result of this section: 

Theorem 2 (Surface code reduction). Let C, C a and Cb be 

surface code systems with the property that there exist qubits 
a, be E(£) such that E(£ a ) = E(£) \ {a} and E(£ 6 ) = E(£) \ 
{b}. A. A a and Ab denote the induced adjacency relations. In 
addition we require 

1. system £ to be stricter than the other systems, i.e. A )p 
A a and A )p Ab- 

2. that there is a leaf-graph G a b, describing an LC- 
equivalent graph state of system £, such that G a b — a 
and Gba — b describe LC-equivalent graph states of C a 
and Cb, respectively. 

3. the surface code systems C a and Cb to be nonlocal. 
Then the surface code system C is nonlocal. 

Proof. In order to show that the surface code system de- 
fined by £ is non-local, one has to find at least one non- 
local edge for each element in [G a t,] LC with respect to the 
adjacency relation A. Let H G [G a 6] LC be an arbitrary 
graph representing an LC-equivalent graph state. Accord- 
ing to Proposition 1, each element in [G a (,] LC has at least 
one subgraph in [G a b — a] LC or [Gt a — o] LC - W.l.o.g. let 
H'a G [G a b — h ]lc ^ e such a subgraph, i.e. H' a C H. Since 
H'a describes an LC-equivalent graph state of system C a (see 
2.) there is at least one edge e = {vi,V2} G E(ff„) C ¥,(H), 
such that the vertices are not adjacent with respect to A Q , 
i.e. e ^ E(Ga q ) (see 3.). Since system £ is stricter than 
C a (see 1.) it follows Ga [E(£ a )] C Ga q and consequently 
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e i E(G A J => e i E(G A [E(£ a )]). Since ui,«2 € E(£ a ) it 
follows furthermore e ^ E(Ga). Thus we found an edge e 
in K(H) that represents a non-local interaction in the larger 
system C. Since H was chosen arbitrarily, this shows the 
non-locality of C. ■ 




Figure 11. Example for the reduction procedure based upon 
Theorem 2. A: The +-shaped pentomino as presented in 
Fig. 2 (A). B-C: By reducing the upper square twice one 
obtains successively smaller systems. D-F: This procedure 
is repeatedly applied to the remaining squares and leads to 
the system depicted in Fig. 2 (C). The "removal" of each 
cdge/qubit is justified by Theorem 2; thus the non-locality of 
the smallest setup (F) implies the non-locality of the whole 
chain of larger setups. 



Using Theorem 2, the +- shaped pentomino (see Fig. 2 (A)) 
can be reduced to the 8-qubit setup depicted in Fig. 2 (C). The 
procedure applied in this particular case is depicted in Fig. 11. 
The peripheral edges are contracted successively until we end 
up with the 8-qubit setup in Fig. 11 (F). The contractions are 
justified by Theorem 2. Each step requires an LC-equivalent 
graph state with leaf-graph G a t as depicted in Fig. 2 (A), 
such that Gab — a describes an LC-equivalent graph state of 
the next (reduced) setup. One can perform these steps easily 
by hand (they are not very exciting!) and, as a consequence, 
verify the non-locality of the -(--shaped pentomino using the 
non-locality of the 8-qubit setup. 



